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Abstract 

The critical exponents and the critical amplitude ratio of the scalar 
model are determined using finite-temperature field theory with aux- 
iliary mass. A new numerical method is developed to solve an evo- 
lution equation. The results are discussed in comparison with values 
obtained from the other methods. 
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1 introduction 



Phase transition is an important phenomenon in particle physics, cosmology, 
and condensed matter physics. Quark Gluon Plasma should be present at 
the heavy ion collision and will give us a lot of valuable information on par- 
ticle physics 0. The investigation into the chiral phase transition suggest 
that a number of flavor may be bounded above [0, |3|]. In cosmology, the 
electro-weak phase transition should be first order for electro-weak baryo- 
genesis [§ § and is investigated attentively @, 0, §, g 0, [TTJ, [T2], |Tj| [TJ]. 
Needless to say, a variety of phase transitions are observed and investigated 
precisely in condensed matter physics. 



Field theoretical approach is essential in order to investigate these phase 
transitions: finite-temperature/chemical-potential field theory [|], [Tj|, pertur- 
bative and non-perturbative renormalization group [[L(| [17], [18|, [19], |20|, 21], ^ 



field theory on lattice |23j, and so on. Temperature can be naturally intro- 
duced by statistical principle using finite-temperature field theory. Not all 
the phase transitions, however, can be investigated by it; the perturbation 
theory, which is the most powerful method at zero temperature, often breaks 
down around the critical temperature because of many interactions in ther- 
mal bath N24L EI| . Indeed the perturbation theory fails, when it is applied to 



either a second-order or a weakly first-order phase transition. 



Drummond et.al.[^6j proposed a new method using an auxiliary mass in 
order to avoid this difficulty We utilized their idea and developed a new 
method to calculate the effective potential. We, then, investigated the phase 
transition of the scalar model using the auxiliary-mass method and showed 
it is second order correctly f27| . It is a great advance in finite-temperature 
field theory, because the phase transition in the scalar model is indicated to 
be first order incorrectly by the perturbation theory with daisy resumma- 
tion ]7|, We note that the method was able to reproduce the result with 



super-daisy approximation 2S 



Since the equation we must solve in the auxiliary-mass method is a non- 
linear partial differential equation for the effective potential, it can not be 
solved analytically and must be solved by a numerical method. It is, how- 
ever, difficult to solve partial differential equations numerically because of the 
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numerical instability JHJ. What is worse, the non-linearity of the equation 
prevents us from using the methods established in case of a linear equation. 
We could not, therefore, make mesh size arbitrary small; The investigation 
in Ref.pTj was not accurate quantitatively. In the present paper, we use an 



improved numerical method given in the appendix, which do not suffer from 
the instability, and get accurate universal quantities. Unlike the rough val- 
ues in p?J , they are beyond the values obtained from Landau approximation. 



The present paper is organized as follows. In the next section we review 
the auxiliary-mass method developed in pTj. In section 13 the effective po- 



tential is shown as temperature varies. We, then, focus on the behaviour 
of it around the critical temperature and calculate the universal quantities. 
These values are compared with values obtained from the other methods. 
Summary and discussion is presented in section [J[ In appendix we devote to 
explain the numerical method we used. 



2 Review of the auxiliary- mass method 



We review the method to calculate an effective potential at temperature 
where the perturbation theory is not reliable [27| . We consider A0 4 theory 
which is defined by the Lagrangian density 




) 2 - -mV 



A 



+ J(f) + C.t., 



where J is an external source function. If m 2 is negative, the scalar field 
4> develops the non- vanishing field expectation value at T = 0. First, the 
effective potential is calculated with a positive mass squared M 2 which is as 
large as the temperature T 2 . This selection of the mass permit us to use the 
perturbation theory without failure, because the loop expansion parameter 
there is XT /M ~ A [0, |3l], [J2|, which is small when the coupling constant A 
is small. Using the perturbation theory, the effective potential is calculated 
as follows, 
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A 74 T 



M 2 72 , '« 74 



drr 2 log 
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(2) 



Here, only the one loop thermal correction is left and the quantum correction 
is neglected, because it should be negligible when the coupling constant A is 
sufficiently small. 

We, then, extrapolate the effective potential @ to the negative mass 
squared[] m 2 = — /i 2 using the following evolution equation, 

dV I79 1 r+i°°+t , r d 3 p 1 

- ■ +ir- d Po ' 



dm 2 2 r 2m J-ioo+e ^ J (2vr) 3 -p 2 + p 2 + m 2 + f <p 2 + IT eP™ - 1 

(3) 

where is an expectation value of the field and II = II (p 2 , —p^, <p, m 2 , r) is 
a full self energy. The thermal correction is left and the quantum correction 
is neglected here, too. Of course, IT can not be calculated exactly; we need 
an appropriate approximation in order to calculate the effective potential 
from (^|). Because the effective potential is a generating function of n-point 
functions with zero external momentum, neglect of momentum dependence 
in II allows us to replace as follows, 

m 2 + ^0 2 + 11(0, 0, 0, m 2 , r) - ~. (4) 

The evolution equation (|3]) can be converted to partial differential equation 
using this replacement as follows, 

~(f> + T~9 / drr 1 ^ = 7 i x • ( 5 ) 




dm 2 T An 2 Jo I WV 

f 2 + w exp 

The effective potential can be calculated by solving the partial differential 
equation (|5[) with the initial condition (|2|). The effective potential has an 
imaginary part below the critical temperature and an analytic continuation 
is done so that this imaginary part is negative pTfl . Since the evolution 



equation (|g) is complicated non-linear partial differential equation, it can be 
solved only by numerical methods. 

1 Hereafter we use unit fx = 1. All dimensionful quantities are measured in the unit. 
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Figure 1: Real part of the effective potential (A = 1). The values of the origin 
are set to zero. A stable point comes to be zero smoothly as temperature 
increases. 

3 Results 



We calculate the effective potential numerically using the method in the Ap- 
pendix. The real part of the effective potential as temperature varies is shown 
in fig.|l|. A stable field expectation value C , where the effective potential has 
its minimum, comes to be zero smoothly as temperature increases. This in- 
dicates that second-order phase transition takes place in this model correctly 



[p3| , p4| . The imaginary part of the effective potential below the critical tem- 
perature is shown in fig.|2|. One can observe a magnitude of it increases as a 
field expectation value decreases; this illustrate that a state with smaller field 
expectation value is less stable below the critical temperature. The critical 
temperature as a function of the coupling constant A is shown in fig.[3[ This 
shows a similar behaviour to the leading result obtained in Ref . , but has 
a slight difference (~ 2%). In the remaining of this section we determine 
some critical exponents: (3,5, 7+/" ,and a. The amplitude ratio x+/X- is 
also determined. The results are summarized in table. [|. 

First, we observe the stable point C carefully. Figure.^ shows C as a 
function of temperature. It decreases monotonically and vanish smoothly as 
temperature increases. We, then, focus on its behaviour near the critical tem- 
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Figure 2: Imaginary part of the effective potential (A = f). The magnitude, 
which shows the instabity of the state, increases as a field expectation value 
decreases. 
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Figure 3: Phase diagram of A0 4 theory. Second-order phase transition is ob- 
served on the boundary. The dots represent values calculated using auxiliary- 
mass method. The dotted line represents the leading result of perturbation 
theory |35 1 . 
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Figure 4: Stable field expectation value as a function of temperature (A = 1). 
It decreases monotonically and vanishes smoothly as temperature increases. 



perature T c and determine /3, which relate a magnetization to temperature 
near T c . This is defined as follows, 

c cx(-rf (r~0,T<T c ) (6) 

where r = (T — T c )/T c . We plot log(0 c ) against log(— r) in fig.||; we fit the 
data to linear function and draw it in fig.|5|. We determine f3 from the gradi- 
ent of it. We find (3 = 0.385. 



Next, we determine the exponent S which is defined as follows, 

0oc JV< = (^)i/* (T = T C ). (7) 

One can derive the following relation from this, 

Foc0 m (T = T C ). (8) 

We show the effective potential at T c in fig.|| We plot log(V) against log(0) 
in fig.0; we fit the data to linear function and draw it in fig.|7|. We determine 
5 from the gradient of it. The result is 5 = 4.0 
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Figure 5: Plot of log(0) — log(— r) (A = 1). The data points are fit to linear 
function. Using its gradient, (3 is determined. 
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Figure 6: Effective potential at the critical temperature (A = f). 
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Figure 7: Plot of log(V) — log(0) (A = 1). The data points are fit to linear 
function. Using its gradient, 5 is determined. 



Then, we determine 7 +// and x+/x— They are defined as follows through 
the susceptibility, 



X ~ dJ 



X = 



dJ 



j=o 



j=o 



X+r 



X-T 



(r~0,T>T c ), 
(t~0,T<T„). 



(9) 
(10) 



To calculate Xi we relate x to the curvature using the following identity 
derived from the definition of the effective potential, 



d(f) 
dJ 



(11) 



We show (jffi) _ as a function of temperature in fig.|[ We also plot 



log(^2") against log(|r|) in fig.0; we fit the data to linear functions and draw 
them in fig.[]. We determine 7 + /~ from the gradient of it and x+/X- from 
the intercepts. We find 7 = 7 + = -y~ = 1.37, x+/X- — 3-4 



8 



d 2 V 

<90 2 




0.001 T 



Figure 8: Curvature at minimum point as temperature varies (A = 1) 



log 



<9 2 V\ 



1e-06 £- 



data points above the critical temperature © 

fit line above the critical temperature --- 

data points below the critical temperature + 
fit line below the critical temperature 




MM) 



Figure 9: Plots of log(|^) - log(|r|) (A = 1). The data points are fit to 
linear functions. Using their gradients, 7+/^ and x+/x~ are determined. 
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Figure 10: Second derivative of the effective potential with respect to tem- 
perature (A = 1). 



Finally, we pay attention to the second derivative of the effective potential 
with respect to temperature, which is proportional to the specific heat C. The 
exponent a is defined as follows^, 



d 2 V 
C oc — — oc r 
or 2 



[r ~ 



0). 



(12) 



This derivative is shown in fig. TD] as a function of temperature. We focus on 
its behaviour around T c in fig. 11 and observe that it blows up there. One of 
the critical exponent a is determined using this. The result is a = 0.12. 

The results are summarized in table .[l| and compared with results obtained 
by various methods. Discussion is presented in the next section. 



2 Though the amplitude ratio of the specific heat can also be defined, it is no determined 
because of the numerical reason. 
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Figure 11: Specific heat C as a function of r around the critical temperature. 
One can observe that it blows up around the critical temperature. One of 
the critical exponent a is determined from this. 



4 Summary and discussion 

The critical exponents and the amplitude ratio were determined using the 
auxiliary-mass method developed in ref. |27| by the improved numerical method 
in the appendix. The results are summarized in table.|I[ We found that A0 4 
theory shows second-order phase transition as it should be. Though the 
critical exponents calculated here do not satisfy the scaling relations, they 
satisfy inequalities of critical exponents. For example, the inequalities given 
by Griffiths |8|, 



7 

7 + (^ + i; 



> W-l), 

> (2-a)(5-l) 



(13) 
(14) 



are satisfied. In the following we compare our result with other's. 

First, the results are compared with the values obtained by the perturba- 
tive finite-temperature field theory with daisy resummation. Since first-order 
phase transition is indicated at one-loop order J7|, the critical exponents 
can not be determined by the perturbation theory. At two-loop order, second- 
order phase transition is observed and the critical exponents are same as those 
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7 


V 





a 


S 


7] 


x+/x- 


F-T 


auxiliary-mass method 


1.37 




0.385 


0.12 


4.0 




3.4 


perturbation 
theory 


1-loop 


* 


* 


* 


* 


* 


* 


* 


2-loop 


1.0 


0.5 


0.5 


0.0 


3.0 


0.0 


2.0 


R-G 


perturbation 
theory [j3t| 


fixed dim. 


1.24 


0.630 


0.325 


0.11 


4.82 


0.317 


4.82 


e-exp. 


1.24 


0.631 


0.327 


0.11 


4.79 


0.349 


4.70 


non 
-perturbative 


sharp cut off IS] 


(1.38) 


0.690 


(0.345) 


(-0.07) 


(5.0) 


0.0 


* 


smooth 
cut off 




(1.32) 


0.660 


(0.33) 


(0.02) 


(5.0) 


0.00 


* 


d' 2 


(1.20) 


0.618 


(0.327) 


(0.146) 


(4.67) 


0.054 


* 


lattice Monte Carlo 31 


1.24 


0.629 


0.324 


0.113 


4.83 


0.027 




experiment ^] 


binary fluids 


1.236 


0.625 


0.325 


0.112 






4.3 


liquid- vapor 


1.24 


0.625 


0.316 


0.107 






5.0 


antiferromagnets 


1.25 


0.64 


0.328 


0.112 






4.9 


Landau approximation 


1.0 


0.5 


0.5 


0.0 


3.0 


0.0 


2.0 



Table 1: Critical exponents and critical amplitude obtained from various 
methods. Since first-order phase transition is indicated, the critical exponents 
can not be determined using finite-temperature field theory (F-T) within one- 
loop order. We note that there are many non-perterbative methods based 
on the renormalization group (R-G) idea which we do not refer here. The 
central values of them are shown. Values in parenthesis are determined using 
scaling relations. 

by the Landau approximation 0. auxiliaryln comparison with these values, 
the results obtained in the present paper are considerably good. 

Second, they are compared with the values obtained by renormalization 
group and by lattice simulation, which agree with greatly. In comparison 
with these accurate values, our results are not very good. These errors are 
probably caused by the replacement ([|). Since this replacement is based on 
the neglect of momentum dependence in II, we have to take into account the 
momentum dependence in order to improve our results [[39]. 



3 We used the two- loop order effective potential calculated in [0. We determined the 
critical exponents from this both numerically and analytically. 
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As mentioned in Sec.[I], the finite-temperature field theory is optimum 
theory in order to investigate phase transitions; it is based on statistical 
principle and can deal with both first-order and second-order phase transi- 
tion. The perturbation theory, however, often breaks down and it prevent 
us from using the finite-temperature field theory. The -mass method enables 
the finite-temperature field theory to be used in various situations. 

We finally express our thanks to T. Inagaki for valuable discussions and 
communications. J.S is supported by JSPS Research Fellowships. 



A Numerical method 



The numerical method, which we use to solve @, is explained in this ap- 
pendix. The partial differential equation (Rf) is written as follows, 



dV l- 2 d 2 V, 

dm? r ; W' 



(15) 



Here, f(x) is the integral in (|5|). First, we make a lattice shown in fig.|12 
The partial differential equation (|5]) is, then, differenced as follows 



30 



Vi 



Vi,j 



Am 2 



Itf + f L(^±M±1 " 2Vw + Vi ~ ld+1 



+(!-«)( 



(A0) 2 

(A0) 2 ■ 



• (16) 



The parameter a decides where the laplacian is evaluated. If a 
selected, the laplacian is evaluated at (a) in fig. 12 



is 

The method with this 



selection is called the explicit method, which we used in [27]. This method is 
simple, because V x j + i is determined only by substituting V x j into the right 
hand side. It, however, suffers from a numerical instability, when smaller 
mesh A0 is chosen [3"0||; therefore, we could not make mesh small in [2~7|. If 
a = 1 is selected, the laplacian is evaluated at (b) in fig.|12|. The method 
with this selection is called the implicit method, which does not suffer from 
the numerical instability at least if f(x) is a linear function |3(| - as far as 
we know, when f(x) is not linear function like our case, not many things are 
known . If a = 1/2 is selected, the laplacian is evaluated at (c) in fig. |12" 
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Figure 12: Lattice used to difference 



The method with this selection is called the Crank-Nicholson method, which 
also does not suffer from the numerical instability at least if f(x) is a linear 
function. What is more, the result converges more rapidly with decreasing 
Am 2 using this method ^Of. Both the implicit and the Crank-Nicholson 
method, however, requires us to solve coupled non-linear equation ([[£]); this 
prevents us from using established method in the case f(x) oc x. 

We developed two methods in order to overcome this difficulty First 
method is based on the Taylor expansion of f(x). The equation fllTf ) is 
rewritten as follows, 



h3 



Am 2 



1,2 f ( Vj+ij - + YizlA 
2 9 ' 1 J { (A0) 2 



Vi+l,j+l ~~ 2 Vij + i + V^_ij + i Vi+lj — 2Vy + 



(A<f>f 



(AW 



(17) 



Since the quantity in the parenthesis behind a is the variation of the laplacian 
per one step, it is small if Am 2 is sufficiently small. We, then, expand f(x) 
around . 

~ 1 ,2 f ( — 2Vij + Vi-ij \ 



Am 2 



+ f 
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+a 



(A0) 2 (A0) 



2 



x f (Yttld — ^ — iJA) -f hiqher order terms (18) 

V ( A </>) 2 / 

This coupled equation is linear with respect to V x j + i and can be solved easily 



The second method is based on an iteration. In order to solve equation 
fllEp, we iterate as follows until a solution is found, 

V n+1 — V- • 1 / V n+1 — ?V n 4- v n 

Am 2 T { J \ 1 (A</>) 2 

+(!"«)( )J- ( 19 ) 

Here, n is the number of the iteration. Note that we can not replace V^" +1 
with Kj+i unlike the Gauss-Seidel method, which is a powerful method if 
f(x) is a linear function |3C|. Next, relaxation method is used in order to 



improve the convergence [30]. Since this procedure is identical with the linear 



case, we only write down the iteration equation without an explanation. 



Am 2 V 2 V ( A< ^) 2 

+(l-a)( 



— 2 ^,j + 



(20) 

Here, the relaxation parameter u; is determined only by experience. The 
results by the two methods agree greatly. In the present paper, the latter 
method is used in order to determine the universal quantities. 
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